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QQ i We study the spectrum of the domain walls interpolating between different chi- 

£\J 1 rally asymmetric vacua in supersymmetric QCD with the SU(N) gauge group 

and including N — 1 pairs of chiral matter multiplets in fundamental and anti- 
fundamental representations. There are always "real walls" interpolating between 
the chirally symmetric and a chirally asymmetric vacua which are BPS saturated. 
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For small enough masses, there are two different "complex" BPS wall solutions 
interpolating between different chirally asymmetric vacua and two types of "wall- 
some sphalerons". At some m = m„, two BPS branches join together and, in some 
interval m* < m < to**, BPS equations have no solutions, but there are solutions 
to the equations of motion describing a non-BPS domain wall and a sphaleron. 
For m > m**, there are no complex wall solutions whatsoever. 



1 Introduction 



f-} • Supersymmetric QCD is the theory involving a gauge vector supermultiplet V 

D , and some number of chiral matter supermultiplcts. The models of this class 

■ attracted attention of theorists since the beginning of the eighties and many 

interesting and non^trivial results concerning their non-perturbative dynamics 
have been obtainedEl. The dynamics depends in an essential way on the gauge 
group, the matter content, the masses of the matter fields and their Yukawa 



■ couplings. 



The most simple in some sense variant of the model is based on the SU(N) 
gauge group and involves N — 1 pairs of chiral matter supermultiplets Si a , 
in the fundamental and anti-fundamental representations of the gauge group 
with a common mass m. The lagrangian of the model reads 
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color and Lorentz indices are suppressed. In this case, the gauge symmetry 
is broken completely and the theory involves a discrete set of vacuum states. 
The presence of N chirally asymmetric states has been known for a long time. 
They are best seen in the weak coupling limit m <C Asqcd where the chirally 
asymmetric states involve large vacuum expectation values of the squark fields 
(sj) S> A-sqcd and the low energy dynamics of the model is described in 
terms of the colorless composite fields Aiij = 2S[Sj. The effective lagrangian 
presents a Wess-Zumino model with the superpotential 

W = 4^£->M (2) 
3 detAl 2 v ' 

The second term in Eq. (^) comes directly from the lagrangian (Q) and the first 
term is generated dynamically by instantons. Assuming Ai^ = X 2 Sij and 
solving the equation dW/dx = (x is the scalar component of the superfield 
X), we find N asymmetric vacua 

/ A 2JV+1\ V2iV 
(X)k= l^QCD ^ k/N ^ Wik/N (3) 

\3 m I 

(the vacua "k" and "k + N" have the same value of the moduli (x 2 )fc and 
are physically equivalent). These vacua are characterized by a finite gluino 
condensate 

(Tr A 2 ) fe = 8^ 2 m( X 2 ) fc (4) 

It was noted recently! that on top of (^) also a chirally symmetric vacuum 
with the zero value of the condensate exists. It cannot be detected in the 
framework of Eq.(^|) which was derived assuming that the scalar v.e.v. and 
the gluino condensate are nonzero and large, but is clearly seen if writing down 
the effective lagrangian due to Taylor, Veneziano, and Yankielowicz (TVY) 
involving also the composite field 

$3 = 32^ ^ < 5 > 
The corresponding superpotential reads 
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Choosing kinetic term in the simplest possible form 
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the corresponding scalar potential is 
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(from now on we set Asqcd = 1)- The potential (||) has N + 1 degenerate 
minima. One of them is chirally symmetric: 4> = x = 0. There are also N 
chirally asymmetric vacua with (x)k given in Eq.(§) and 
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The presence of different degenerate physical vacua in the theory implies 
the existence of domain walls — static field configurations depending only 
on one spatial coordinate (z) which interpolate between one of the vacua at 
z = — oo and another one at z = oo and minimizing the energy functional. As 
was shown inu, in many cases the energy density of these walls can be found 
exactly due to the fact that the walls present the BPS-saturated states. 

The energy density of a BPS-saturated wall in SQCD with SU(N) gauge 
group satisfies a relationB 
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where the subscript ±oo marks the values of the gluino condensate at spatial 
infinities. 

Bearing Eqs. (O, 0) in mind , the energy densities of the BPS walls are 
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for the real walls interpolating between a chirally asymmetric and the chirally 
symmetric vacua (we call them "real" because the corresponding gauge field 
configurations are essentially real) and 
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for the complex walls interpolating between different chirally asymmetric vacua. 
The RHS of Eqs.(|lC|-|l2|) presents an absolute lower bound for the energy of 
any field configuration interpolating between different vacua. 

The relation ( |I(]| ) is valid assuming that the wall exists and is BPS- 
saturated. However, whether such a BPS-saturated domain wall exists or 
not is a non-trivial dynamic question which can be answered only in a specific 
study of a -particular theory in interest. This talk presents a brief review of 
our papersEra where such a study has been performed. 

The results concerning the BPS walls were obtained by solving numerically 
the first order BPS equations 



d z cj) = e i5 dW/d<t>, d zX = e w dW/dx 
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associated with the TVY lagrangian. The value of 8 to be chosen in Eq.([L3|) 
depends on the wall solution we are going to find. To fix it, note that the 
equations (13) admit an integral of motion: 
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Indeed, we have 
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The real wall connects the vacua with W = and W = — N/2(Ai 



These boundary conditions are consistent with Eq.(|f4|) only if e lS = ±1 (the 
sign depends on whether the walls goes from the symmetric vacuum to the 
asymmetric one when z goes from —oo to +oo or the other way round). For 
the complex walls, the boundary conditions and the conservation law (14) 



imply 5= f ±f. 

To study the spectrum of the domain walls which are not BPS-saturated, 
one has to solve the equations of motion which are of the second order, and, 
technically, the problem is a little bit more involved. 



2 Real Walls 

The BPS equations (|T^) with the superpotential (|J) and the kinetic term (Jfy 
have the form 
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Figure 1: Real walls (N = 3). The parametric plots for m = 1 (solid line), m = .1 (dashed 
line), and m = 10 (dotted line). <fi an d X are measured in units of R t and p», respectively. 

(O' = d z O). To find the wall interpolating between <jf> = x = 0at2 = — oo 
and (/) — R^,x — P* at z = oo, we have to choose 5 = 7r (or <5 = for the wall 
going in the opposite direction) . With this choice and the boundary conditions 
given, the solutions 4>{z) and x( z ) are going to be real so that we have a simple 
system of just two first-order differential equations. 

The dynamics of this system is essentially the same for any N. The solution 
exists for all masses. For large m, the heavy matter field can be integrated 
out, and we arrive at the BPS equation for the supersymmetric giuodynamics 



The solution of this equation with the boundary conditions </>(— oo) = 0, </>(oo) = 
R* can be expressed into integral logarithms. 

An analytic solution of the system ( |l5| ) can be found also for small masses. 
In the intermediate range of masses, the solution has to be found numerically. 
Parametric plots in the (</>, \) plane for N — 3 and different values of m are 
drawn in Fig. ^. 

3 Domain Walls in Higgs Phase 

The lagrangian (||) can be obtained from the full TVY lagrangian (||) in the 
limit of small masses if integrating out the heavy gauge under the assumption 
that the Higgs expectation value of the scalar matter fields is large. In this 
approach, chirally symmetric vacuum is not seen and there are only complex 
walls which turn out to be BPS saturated (we have seen, however, that the 
chirally symmetric vacuum and the real walls associated with it exist actually 
for all masses no matter how small they are). To find the profile of complex 
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walls in the Higgs phase, we have got to solve the BPS equations 

X ' = e --(iv-2)/2iV^ 
dx 

where W is the superpotential @), with the boundary conditions 

x(-oo) = p*; X (oo) - p^l N - (18) 



It is convenient to introduce polar variables \ — pe la . The equations (17) 
acquire the form 

P' = (iV-l)|mpsin 7 - sin[ 7 (iV - 1)] 



V = 2lV-1)^/;m-os-. :- JLjcosMN-l)^ 



where 7 = 2a — 7r/iV changes from 7 = —n/N at z = —00 to 7 = 7r/7V at 
2 = 00. For N = 2, the solution is analytic: 

p(z) = 

tan 7 (z) = sinh[4m(z - z )] (20) 

or in the complex form: 

1 + j e 4m ( z - z o) 

X (z) = p. - — - (21) 
V 1 + e 8m ( z_2; o) 

(zo is the position of the wall center). For TV > 3, the solutions can be found 
numerically. The profiles for the ratio r{z) — p{z)/p Sf in the interval z$ = 
< z < 00 [it is a half of the wall, another half being restored by symmetry 
considerations: p(— z) = p(z)\ with different N = 3, 5, 10 are presented in Fig. 

I 

We see that the dependence p(z) is not flat anymore but displays a bump 
in the middle. To understand it, remind that the system ( fl9| ) has the integral 
of motion (|l4|). In our case, it amounts to 

m(N-l) 2 ^ 2 AT (Am N - 1 \ 1,N 



V cos 7 - ^-^^ cos[ 7 (iV - 1)] = _ j cos - (22) 

as follows from the boundary conditions ( )l8|) w ith S = n/N — n/2. In the 
middle of the wall, 7 = 0, and the condition (|22|) implies 

^-^-^T) = Nco S ^ (23) 
6 




Figure 2: BPS walls in Higgs phase for N = 3 (dashed line), N = 5 (dotted line) and 
N = 10 (solid line). 



( x = r(0)). It is not difficult to observe that the real root of the algebraic 
equation (^) is slightly greater than 1 for N > 3. When N is large, x — 1 
tends to zero <xl/N. 

4 Two BPS solutions and Phase Transition in Mass 

When m ^ 0, the gauge degrees of freedom associated with the superfield <& 
do not decouple completely and should be taken into account. The full system 
([La) of the BPS equations for the complex domain walls has the form 



P' = (N-l) 
a' = (N-l) 



mp8ia(2a - n/N) - ^ sin(3/3 - n/N) 



||- 008(3^9 - n/N) 



m cos(2a — n/N) 

R' = -2R 2 [sin(3/3- n/N)bx{R 3 p^ N ~^) + cos(3(3 -n/N)[3(3 + 2a(N - 1)]J 
/?' = 2i?[-cos(3/3-7r/iV)ln(i? 3 / 9 2(Ar_1) )+sin(3/3-7r/iV)[3/3 + 2a(iV-l)]] 

where, as in the previous section, we have chosen S = it /N — n/2 and intro- 
duced the polar variables x — P eta > 4> — Re 113 ■ One should solve the system 
(pif) with the boundary conditions 
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R{oo) -- 
2{N - 



p(-oo) = p(oo) = p*; i?(-oo) 

a(-oo) = /3(-oo) = 0; a(oo) = n/N; /3(oo) 

In the limit m — > 0, we can just freeze the heavy variables 
13 = -2a(N - l)/3; R = p^" 1 )/ 3 
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in which case the first two equations in Eq.(|24|) reproduce the system (19) 
studied in the previous section. When mass is nonzero, but small enough, it is 



possible to develop a systematic Born-Oppcnhcimcr expansion and to find the 
profile of the wall as a series over the small parameter m (or rather m( 2JV+1 ^ 3JV 
in our case). Not going in details 0, we just present here for illlustration how 
R(0), the absolute value of the field 4> in the center of the wall, is changed with 
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where R (0) oc m ( Ar ~ 1 )/( 2JV+1 ) is related to p (0) = xp* as in Eq.@, and 
Pq(0) was found earlier when we studied the solutions of the reduced system 
in the Higgs phase. 

When to is neither too large nor too small, the Born-Oppcnhcimer ap- 
proximation does not apply and we are in a position to solve the full system of 
4 equations (|24| ) numerically. We did it for N = 2, 3, 4. The results turned out 
to be rather nontrivial and surprising. Namely, there turned out to be not one 
but two different BPS solutions. But they exist only for small enough masses; 
if the mass exceeds some critical value to*, the BPS equation system has no 
solution at all. A kind of phase transition in mass occurs ! The values of the 
critical mass m* have been found numerically: 



m 



SU(2) 



4.67059 . 
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.28604. 
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0.07539... (28) 



In Fig. ||, we plotted the dependence of r\ = R(0) /Rq(0) on m for both branches 
. We see that, at to — m*, two branches are joined together. This is the reason 
why no solution exists at larger masses. 

The upper BPS branch approaches the Higgs phase solution for the small 
masses. Speaking of the low branch, it approaches two widely separated real 
BPS walls for N = 2 passing through the chirally symmetric vacuum in the 
middle, but for N > 3 it presents a new nontrivial solution with the field 
values in the middle of the wall passing at some finite distance from the chirally 
symmetric vacuum even in the limit to — > 00. 

As is seen from Eq. (p8|), the value of m* expressed in the units of 
Asqcd drops rapidly with N. One can understand it looking at the Born 
Oppenheimer result ( p7| ) Indeed, the critical mass to* and the value of the mass 
where the Born-Oppenheimer expansion in Eq.(^7|) breaks down should be of 
the same order. For large N, the expansion parameter is n ~ A 2 (m/i?*) ~ 
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Figure 3: The ratio r) = -R(O) /Ro (fl) as a function of mass for the SU (3) theory. The dashed 
line describes the analytic result ( [27] ) valid for small masses. 

iV 2 m 2 / 3 . Assuming k ~ 1, we arrive at the conclusion that m* falls down as 

m*(iV) cx TV" 3 

in the limit iV — > oo. 

5 Solving equations of motion 

Every solution of the BPS equations ( |l3| ) is also a solution to the equations of 
motion, but not the other way round. It turns out that nontrivial non-BPS 
complex wall solutions exist. To find them, we have to solve the equations of 
motion for the potential (^) with the kinetic term \d(f>\ 2 + |<9x| 2 - As earlier, it 
is convenient to introduce the polar variables x = P e%a '■> — after which 
the equations of motion acquire the form 

) 



(29) 
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Figure 4: The ratio r) = R(0) / Rq(Q) for the solutions of the equations of motion as a 
function of mass for the SU (2) theory. The solid lines describe the BPS solutions and the 
dashed lines describe the non-BPS wall and the sphalerons. 



where L = bjjii 3 /? 2 ^ -1 )], 0+ = 3j3 + 2(N- l)a, /3_=3/3-2a. The boundary 
conditions (|25|) are assumed. 

The phase space of the system ( p9[ ) is 8-dimensional and the boundary 
problem (which we solved by shooting method) is somewhat more tricky than 
for the BPS case. In particular, the trajectories become more and more unsta- 
ble as N grows. Finding the trajectory going from one vacuum to another is 
something like lancing the ball to roll along a razorblade mountain ridge from 
one top to another. Our computer managed to do it for N = 2 and N = 3. 

All the solutions obtained are presented in Fig. ^ (N — 2) and in Fig. 
H (N = 3) where the parameter R(0) is plotted as a function of mass [we 
normalized R(0) at its value R o {0) = R* (TV = 2) and R (0) = .918... it, 
(N = 3) for the upper BPS branch in the limit m — ► ]. For small masses, 
there are several solutions. We obtain first of all the BPS solutions studied 
before (solid lines in Figs. ||, ||). We find also two new solution branches drawn 
with the dashed lines. We see that, similarly to the BPS branches, two new 
dashed branches fuse together at some value m = m*, > to* . No solution for 
the system ( p9f ) exist at to > m**. 

We see that the pattern of the solutions for N = 2 and N — 3 is qualita- 
tively similar. To understand it, assume first that to < to* and draw the energy 
functional E for field configurations with wall boundary conditions minimized 
over all parameters except the value of R(0) which is kept fixed (see Fig. |^a 
). For very small i?(0), our configuration nearly passes the chirally symmetric 
minimum and the minimum of the energy corresponds to two widely separated 
real walls. Thus E[R(0) = 0] = 2e r with e r given in Eq.([ii"|). Two minima in 
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m 

Figure 5: The same for N = 3. 

Fig. |£|a correspond to BPS solutions with the energy e c . They are separated 
by an energy barrier (for illustrative purposes, the hight of the barriers is ex- 
agerated). The top of this barrier (actually, this is a saddle point with only 
one unstable mode corresponding to R(0), in other words — a sphaleron) is a 
solution described by the upper dashed line in Figs. [|, ||. The lower dashed 
line corresponds to the local maximum on the energy barrier separating the 
lower BPS branch and the configuration of two distant real walls at R(0) = 0. 

At to = to*, two BPS minima fuse together and the energy barrier sepa- 
rating them disappears. The upper sphaleron branch coincides with the BPS 
solution at this point. When m is increased above m», the former BPS mini- 
mum is still a minimum of the energy functional, but its energy is now above 
the BPS bound (see Fig]|b). The corresponding solution is described by the 
analytic continuation of the upper sphaleron branch. The lower dashed branch 
in the region to* < to < m** is still a sphaleron. At the second critical point 
to = m**, the picture is changed again (see FigJ|c). The local maximum and 
the local minimum fuse together and the only one remaining stationary point 
does not correspond to an extremum of the energy functional anymore. At 
larger masses, no non-trivial stationary points are left. 

There is one distinction, however. For N = 2, the second critical mass 

For N = 3, in contrast to the case N = 2, the existence of such a barrier could not be 
established from the BPS spectrum alone. The matter is that, while 2e r = e c for N = 2 and 
the presence of the maximum is guaranteed by the Roll theorem, 2e r > e c for N = 3, and 
one could in principle imagine a situation where E falls down monotonically when R(0) is 
increased from zero up to its value at the lower BPS branch. Our numerical results strongly 
suggest, however, that the energy barrier ( though a tiny one) is present. 
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Figure 6: Illustrative profiles of the energy functional vs. R(0). 



to** ~ 4.83 is very close to m* and the energy of the non-BPS wall at to — m** 
deviates from the BPS limit by less than 10~ 5 ! When N — 3, to** « 3.704 
differs essentially from to* and the energy of the wall at to = to**, being very 
close to the energy of two widely separated real walls 2e r , differs from the BPS 
limit e c for the complex walls. One can speculate that, while to* drops rapidly 
as N grows, the value of to** is roughly constant in the large N limit. To be 
precise, we did not perform a careful numerical study for N > 3 and cannot 
therefore exclude the possibility that, for large enough N, the second critical 
mass disappears altogether and the non-BPS complex wall solutions exist for 
any mass. Our bet, however, is that it does not happen. 

6 Discussion. 

Our main result is that, while the real BPS domain walls connecting the chirally 
symmetric and a chirally asymmetric vacua are present at all masses, the 
complex BPS walls interpolating between different asymmetric vacua exist 
only for small enough masses to < to*, m* being given in Eq. (p8|) . A kind of 
phase transition associated with the restructuring of the wall spectrum occurs. 



Furthermore, also a second critical mass m** exists beyond which no com- 
plex wall solution can be found whatsoever. That means in particular that 

c Nccdless to say, it is not a phase transition of a habitual thermodynamic variety. In par- 
ticular, the vacuum energy is zero both below and above the phase transition point — 
supersymmetry is never broken here. Hence E vac (m) = is not singular at m = m*. 
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no domain walls connecting different chiraly asymmetric vacua are left in the 
pure supersymmetric Yang-Mills theory corresponding to the limit m — > oo, 
and only the real domain walls connecting the chirally symmetric and a chirally 
asymmetric vacuum states survive in this limit. Note that that contradicts a 
recent assumption by Witten lI that it is complex rather than real domain 
walls which are present in the pure SYM theory (Witten discussed them in the 
context of brane dynamics). 

One has to make a reservation here: our result was obtained in the frame- 
work of the TVY effective lagrangian (^) whose status [in contrast to that of 
the lagrangian (||)] is not absolutely clear: the field $ describes heavy degrees 
of freedom (viz. a scalar glueball and its superparnter) which are not nicely 
separated from all the rest. However, the form of the superpotential (g) and 
hence the form of the lagrangian for static field configurations is rigidly dic- 
tated by symmetry considerations; the uncertainty involves only kinetic terms. 
It is reasonable to assume that, as far as the vacuum structure of the theory is 
concerned (but not e.g. the excitation spectrum — see Reftl for detailed dis- 
cussion), the effective TVY potential (||) can be trusted. A recent argument 
against using Eq. ^ that the chirally symmetric phase whose existence follows 
from the TVY lagrangian does not fulfill certain discrete anomaly matching 
conditions E2l is probably not sufficient. First, it assumes that the excitation 
spectrum in the symmetric phase is the same as it appears in the TVY la- 
grangian which is not justified. Second, it was argued recently that the TVY 
lagrangian describes actually all the relevant symmetries of the underlying the- 
pry and the absence of the anomaly matching is in a sense an "optical illusion" 

It would be very interesting to repeat the numerical study performed here 
for the effective lagrangian with the same TVY superpotential but a dif- 
ferent kinetic term. Probably, it makes sense to start with the VY lagrangian 
for the pure SYM theory and try tentatively for example 



(One cannot write instead a general Kahler term J d 4 8F(<&,$>). It would in- 
volve a dimensionful scale and would spoil the anomaly matching condition 
d£ c g/d(\nA) ~ G 2 .). Our guess is that, for the VY lagrangian with the ki- 
netic term (pQ), the real walls talking to the chirally symmetric vacuum would 
also be present and the complex walls would be also absent. For a generalized 
TVY lagrangian, the particular values of m* and rn*» would be different, but 
the qualitative pictures in Figs. 0, F| would not change. 




(30) 
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